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We discuss how to implement the source terms in Quantum Chromodynamics (QCD) respecting 
gauge invariance and noncommutativity of color charge density operators. We start with decom¬ 
posing the generating functional of QCD into constituents projected to have specified color charge 
density. We demonstrate that such a projection leads to the gauge invariant source terms consisting 
of a naive form accompanied by the density of states which cancels the gauge dependence. We 
then illustrate that this form is equivalently rewritten into a manifestly gauge invariant expression 
in terms of the Wilson line. We confirm that noncommutativity of color charge density operators 
is fulfilled in both representations of the source terms. We point out that our results are useful 
particularly to consider the problems of the quantum evolution in high energy QCD. 

PACS numbers: 12.38.-t,12.38.Aw 


I. INTRODUCTION 

The first-principle calculations based on the fundamental theory of the strong interaction, i.e. Quantum Chromo¬ 
dynamics (QCD), are hindered by at least two major obstacles. One difficulty originates from the fact that the theory 
is nonlinear in gauge fields and it is hard to find a solution even at the classical level |l|. Another arises from the 
nonperturbative nature, which indicates various differing aspects of QCD depending on the context. The most obvious 
manifestation of the nonperturbative nature is, of course, that the strong coupling constant g is substantially large. 
The strong coupling constant, however, runs with the relevant momenta conveyed by gluons so that at high enough 
energies the perturbative QCD (pQCD) is expected to be reliable owing to the asymptotic freedom. 

Even when the perturbation theory seems plausible, nonlinearity must be properly taken into account once the 
parton distribution becomes large, i.e. the saturation effects are important in the small-x region @J. Bjorken’s x is 
the longitudinal momentum fraction x = p + /P + in the infinite momentum frame where the target hadron moves at 
P+ = oo. Recently powerful tools to look into small-x wee partons have been developed from the deliberation of this 
saturation picture, which leads us to an effective theory of small-x partons described by the classical equations of 
motion. The classical description is validated because dense partons (or gluons A£) at small x are expected to behave 
like the classical fields created by the source p a that are brought in by the partons with larger x nan In the 
classical model the sources are just regarded as the color charge density carried by valence quarks inside the target 
hadron provided that the target is a large nucleus. The classical gluon fields as strong as A 1 ', ~ 1/g in this picture are 
called the color glass condensates. 

The classical model is the lowest-order approximation of QCD at small x given large p a stemming from larger-x 
partons. One can improve the approximation by taking account of quantum corrections around A%. The theoretical 
framework at the one-loop level has been well organized and the central result is now known as the Jalilian-Marian- 
Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) equation i a s a m m 0 The JIMWLK equation is 
widely accepted, and yet, the treatment of the source terms still leaves some subtleties, which we are addressing here. 

First of all, the gauge invariant generalization of the source terms is not unique. The minimal requirements are 
that the source terms are gauge invariant and are reduced to the naive (gauge variant) form ~ tv[pA~\ not to affect 
the classical equations of motion. The simplest choice would be either ~ tr[pW[A - ]] as adopted in Ref. [d] or 
~ tr[pln W[A - ]] as proposed in Ref. 13] with anticipation from Wong’s equations, where W[A - ] is the Wilson line 
in terms of A~ in the light-cone temporal direction. One can also consider as many variants as one likes which satisfy 
the minimal requirements. So far, concerning the quantum evolution at the one-loop level, all lead to the equivalent 
results, but in principle, the choice should be uniquely determined according to the approximation adopted. We will 
see later that the eikonal approximation for larger-x partons naturally gives rise to the source terms in the form 
~ tr[pln JF[A~]] which has some intriguing properties that tr[pW[A~]] does not have. 

Second of all, the color charge is not commutative at the operator level, that means a combination of different 
colors can make another color. In the classical model the nontrivial commutation between color charge density is not 
considered since p a is assumed to be large enough in the dense regime and the noncommutativity can be ignored. Such 
an approximation is no longer valid in the dilute reg ime, however. Thus it has been proposed that the introduction of 
the Wess-Zumino term can handle this issue [lj,[l5| whose necessity is also understandable from the gauge invariance 
of the source terms jlfil . In later discussions we will clarify that under the eikonal approximation it is the gauge 
variant density of states which plays an equivalent role as the Wess-Zumino term in the sense as discussed in Ref. fl6l . 
Also we will show that the source terms expressed by IT[A~] have already taken care of noncommutativity properly. 
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II. SOURCE TERMS 

The QCD generating functional can be regarded as a sum of all the contributions with distinct color charge density 
p a (x). In this sense we can regard the generating functional as the partition function in grand canonical ensemble 
with zero chemical potential. We will here manage this using the projection operator which imposes a constraint onto 
the functional integral. We essentially extend the formulation of the canonical description of QCD with respect to 
the quark number 0 , 0 to the case of color charge density. 

We denote by Q the Gauss operator which projects on states with a given color charge density p and start with 
decomposing the generating functional as 

Z = z CE[p}= [vAVB g{A,B;p\e iSYM[A]+iSmM ° r[A ’ B] , (1) 

where the light-cone gauge is implicitly assumed for the gauge fixing prescription. We could call Zqe [p\ the partition 
function in the canonical ensemble in the same sense as in Ref. |18| . namely, the partition function for states with 
a given distribution of the color charge density. In our notation B represents the fast-moving matter field bearing 
larger-x which generates the source p a (x) for small-x partons. In the classical model of the color glass condensate 
which is relevant at not very small x it is the valence quark field ip a that B represents, and at small x where quantum 
evolution is crucial B is the fast-moving gluon field A%. We will argue respective cases in order below. 

From the decomposition CD we can define the source action Sw after integrating over B as follows; 

*= £ 

fp(x)} 

All we have to do to identify S-w[A\ p] are thus to write down Q\A , B ; p] and to integrate B out explicitly. 


T>A e iS ™i\A\+iSw[A\p\' 


( 2 ) 


A. Fundamental Representation 

Let us first consider the case that B is the valence quark field. Here we can write the projection constraint explicitly 
as follows; 


Sf[p\ = n^ a ( x ) +#7 + ^W] = fv<t> exp 

a x 


d A x(f) a {x){p a {x) + g-i/j^ + t a ijj{x)} , 


(3) 


where t a, s are the SU(3) algebra in the fundamental representation which satisfy the commutation relation [t a ,t b ] = 
ifabcfc an j are normalized as tr t a t b = 5 ab . The integration measure associated with the (/^integration is invariant 
under color rotations by definition. In the context of the canonical ensemble <p a is often interpreted as the imaginary 
chemical potential. The projection operator actually acts as the Laplace (Fourier) transformation from the (imaginary) 
chemical potential to the density p a . 

The fermionic fields are now supposed to obey the antiperiodic boundary condition with a period 2T, which is 
inspired from knowledge on the Unite-temperature field theory and is known to be suitable for describing the fermionic 
particle distribution. We shall pick up the quark propagation only in the light-cone temporal x + direction as quarks 
are fast moving. This manipulation corresponds to the eikonal approximation since neglecting spatial derivatives di 
means that we drop any hopping or recoil in the spatial directions. The action for matter is thus given by 


^matter [V’,'0] = J d*X 4>i'y + (8 - igA a t a )lp . (4) 

If we use the notation as in the derivation of the JIMWLK equation 00 , A in the above expression would be 
the soft gluon 5A~ actually. There must be its dual , d + — iga a t a , in the presence of the classical background gauge 
field Aa = a a in the covariant gauge if we consider not only the target but also the projectile hadron. We shall leave 
discussion on the projectile to future work and neglect the latter throughout this paper. 

Then the structure of the Dirac matrices becomes trivial because only 7 "*" appears and the integration over the 
quark fields in the presence of the projection insertion results in the Dirac determinant, 


My[A ,<j>\ = TT det 7 + Det d — ig(A a + </> a )f 

A A Dirac 


( 5 ) 
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Here Det stands for taking the determinant in time as well as in the color indices. It should be noted that 4> a always 
appears at the same place as A~. Therefore we shift the integration variable in the following way, 

4>a -» (fa - A~ , (6) 


in order that the determinant becomes a function of (f> a alone. 

It is quite important to note that Det[<9 _ — ig(f> a t a ] = Det[U($ _ — ig<f> a t a )V'} indicates that A4 f[</>] is invariant 
under the gauge transformation of <j> a , 

<t> a r^v(4at a -^vi. ( 7 ) 

The determinant can be evaluated most easily if V is chosen such that (f a t a is static (independent of x + ) and diagonal 
in color (3x3 fundamental) space. Once we eventually get an expression invariant under the gauge transformation 
of (f> a , it can be rotated back to general (j) a - We see that a particular choice of the periodic form, 

VHx) = W* +t _ T (x) [W^_ T (x)\ (x++T)/2T S\x ), (8) 


suffices for our purpose. Here we defined the Wilson line at each spatial point x = (x ,x t ) written as 


W^_ T (x) = V exp 



(9) 


and S^(x) is a color matrix which diagonalizes the zero-mode. It should be noted that V must be periodic, otherwise 
the antiperiodic boundary condition for if is broken. In © the second part [WP_ T ]( X + T i/ 2T is placed to restore the 
periodicity that is violated by the far left part W% + _ T . Then, substituting © into 0, one acquires the transformed 
<t> a t a as 


(j) a t a -> diag(ai,a 2 ,a 3 ) = —^-5(ln , (10) 

i2Tg 

where W^ _ T . In general In is not diagonal in color space and S i is chosen so as to rotate In to a 

diagonal matrix diag(ai, a 2 , a 3 ) with a constraint a 3 + a 2 + a 3 = 0. We do not have to know an explicit form of S'. 

Under the antiperiodic boundary condition in x + with the period 2T, the determinant is evaluated on the basis 
e lUnX+ with the discrete (Matsubara) frequency uj n = (2n+l)n/2T. Then by means of the formula, 11)7=0 [1 — a; 2 /(2n + 
l) 2 ] = cos(7ra;/2), we can easily take the determinant with respect to the x + direction to have 

M k [0] = TT det [d~ — igai ] 4 det [c> - — *ga 2 ] 4 det — iga ^\ 4 

x 

= [cos(Tgai)cos(T3a 2 )cos(T3a 3 )] 4 

X 

= J|[2 + trW^(f)+trIU^(f)] 4 , (11) 

X 

where in the last equality we used a\ + a 2 + a 3 = 0 and 

e z2Tg ai + e z2Tga 2 + T g a 3 = tr ex p [i2Tg diag(ai, a 2 , o 3 )] = tr exp [5(ln W^S^] = trW^ (12) 

to make it take an obviously gauge invariant form tl9ll . The physical meaning is now transparent; no phase factor 
with the weight 2 corresponding to the propagation of no quark and color-singlet three quarks, W^ corresponding to 
the propagation of one quark, and W^ the propagation of two quarks which is equivalent with that of one antiquark. 
The power 4 comes from the determinant in the Dirac indices, namely, spin and quark-antiquark degeneracy. We 
would draw attention to the point that A4 f is written only in terms of trW^, which guarantees the color commutation 
relation as we will see in Sec. m 

We shall present an alternative expression of CD for later convenience. The product over x can be expanded to result 
in the sum over all the charge distributions. That is, ignoring irrelevant overall constants and using a 3 + a 2 + a 3 = 0, 
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we can write 




J ] [1 + e ~ l2Taa '] 4 [1 + e ~ i2T9a2 ] 4 [1 + e ~ t2T9a3 ] 4 

X 


^ ' 4^'n 1 (x) 4Cn 2 (x) 4C ns ^ 

x m (x) ,U 2 (x) ,n3 (x)=0 


i2Tg[n 1 (x)a 1 (x)+ri2(x)a2(x)+n3(x)a3(x)] 


J2 >V F [p] e-f tr ^ sin ^ st ] 

{P(5T)} 


(13) 


where 4 C n = 4!/(4 — n)\/n\ which is the combinatorial weight defining the weight function Wf- From the second to 
the third line we identified p = diag(|gni, \gn. 2 , \gn 3 ) and used 1101) . The explicit form of Wf has been discussed 
in the random walk picture [2Cj, from which not only the Gaussian form in the classical model but also the C-odd 
terms relevant to Odderon arise j2lj . The information of the target size should enter Wf through the measure which 
is necessary to convert the product over x to the exponential of the integral over x. 

In 1131) the summation over p(x ) seems to be taken only with respect to the diagonal components instead of eight 
color charges. This is because (f> a t a was chosen to have only the diagonal components and thus the summation with 
respect to off-diagonal components of p would, if any, result in irrelevant constants. Actually, even though it merely 
multiplies irrelevant constants when (j> a t a is diagonal, one should consider that contains the summation with 
respect to off-diagonal components as well. This is required from the property, Mf[^ a t a ] = Mf\U cpa^U^}, which 
follows for any static U(x) (see 0). Such a transformation in m causes the color charge rotation as p —> U^pU. 
This invariance can be restored if the summation is symmetrically taken over all components of p and moreover Wf 
is augmented to have rotational symmetry in color space, Wf [p] = Wp[UpW], which enables us to eliminate S(x) 
appearing in m and m by redefining S* 1 pS —> p. After all we reach the general form, 

Mf[(!>} = ^ Wf[p] expj- [d 3 xti-[p(x)lnW‘t > (x)]\, (14) 


where {p{x)} distributes not only over the diagonal but the off-diagonal components also. 


B. Adjoint Representation 


In the case when fast-moving partons are gluons, the projection constraint takes a form of 


G,iP\=ims[P a (*) + D z b K + ] 

a x 


J V(j> exp 


d i xc!> a {x){p a {x) + DfF^} 


(15) 


Then (f> a is always accompanied by A a again and we first perform the shift as in the previous case. 

The full integration over the gluon fields is not feasible because of the presence of nonlinear higher-order terms. We 
shall thus perform the one-loop calculation, that is, we divide the gluon fields into fast gluons a% with large p + and 
slow gluons A 9 with small p + , and expand the Yang-Mills action in terms of a 9 up to the quadratic order, and then 
perform the Gaussian integration with respect to a 9 . 

Under the eikonal approximation that p + is much larger than other scales we shall consider in the gluon propagator 
only the terms involving d + D~ = d + (d~ —ig(j> a T a ) with the SU(3) algebra T a in the adjoint representation. We make 
use of the freedom similar to 0 with adjoint matrices to make 4> a T a static and aligned to the Cartan subalgebra 
(chosen as T 3 and T 8 which are not diagonal in the adjoint representation) as 


4>aT a -> foT 3 + </> 8 T 8 


0 —ai2 0 0 0 0 0 O' 

o 12 0 0 0 0 0 0 0 

0 000 00 00 

0 0 0 0 o 3 i 0 0 0 

0 0 0 -031 0 0 0 0’ 

0 0 0 0 0 0 -o 23 0 

0 0 0 0 0 o 23 0 0 

0 000 00 00 


(16) 


where we defined Oi, 02 , and 03 = — 01 — 02 as <(>3 = 01 — 02 and cf> g = \/3(ai + 02 ) in accord to m- We further 
defined 012 = Oi — 02 , 023 = 02 — 03 , and 031 = 03 — oi. We can regard 012 as the charged gluons in the color 1-2 
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sector, 023 i n the 6-7 sector, and 031 in the 4-5 sector likewise. This matrix structure is general also for the SU (N c ) 
case, which would be clearer in the ladder basis in color space [22l[23|. 

After some algebra we find that the integration over the gauge field yields the — § x 2 = —1 power of the determinant 
of the kinetic part d + (d~ — ig<j) a T a ) in the light-cone gauge, A + = 0, where 2 comes from the number of transverse 
gluons. Interestingly, Det(9“ — ig<t> a T a ) is identical with the Haar measure (or the Faddeev-Popov determinant 0 ), 
which eventually leads to the Vandermonde determinant, 


M~ A 1 [(t>\= n Det (d~ 

X 


ig<t> a T a )= nn 

x i<j 


e i2Tga.i _ £ i2Tgaj 


sin 2 ( Tgai 2 ) sin 2 (T ga 23 ) sin 2 ( Tga 31 ) 

X 


(17) 


under the periodic boundary condition with the period 2 T. We note that the periodic boundary condition is appro¬ 
priate to realize the distribution of bosonic particles. One can readily confirm CE3 in the same way as in EJ by 
using the discrete frequency tu n = 2mr/2T for bosons and the formula, n^=i[l — x 2 /n 2 } = sui(ttx)/ttx. We dropped 
irrelevant overall factors. Here the zero-mode of the 3-3 and 8-8 components of m in the frequency product seems 
to render the whole determinant vanishing at a first glance, but it should be replaced by an irrelevant nonvanishing 
constant which in fact comes from the neglected spatial derivatives ~ d 2 . From another point of view, one can say 
that the singular zero-mode should be removed by the complete gauge fixing on the transverse zero-modes 1251 . In 
any case what we need to know is the dependence on (j> a T a (or at) and the zero-mode that has no effect on such 
dependence is safely dropped. 

It is possible to express m by the traced Wilson line in the fundamental representation in the same way as the 
previous case of the fundamental representation, that is [llil ] 


Ma[4>}= J] 


1 - ^|trW^| 2 - A^trW^) 3 - ^_|trW*| 4 


(18) 


It is also instructive to reexpress the results in a form similar to HJ that is convenient for later discussions. By 
using sin 2 (Tga) = §(1 — cos(2 Tga)) and expanding its inverse in terms of cos(2 Tga) one may write the determinant 
(not the inverse) as 


OO 

mai-h =n e [5 

L2 —0 
00 

E[* 


e i2Tgai2 e ~i2T ga\ 2 ^ 


x ni2=0 

X 


1 / e i2Tga23 _|_ g— i2Tga 2 3 s ^ 


77-23—0 


>] E [y 


e i2Tga.3i _|_ g — i2Tga.3\ ' 


n 31 =0 


(19) 


The intuitive meaning is now clear; there are 77 , 12 , ^ 123 , and 7731 gluons of 0 , 12 , 023 , and 031 respectively at each x 
and each takes either positive or negative charge with equal probability. The product over x can be written with an 
appropriate weight function Wa in the form of summation over the charge distributions like in O- If the product 
is expanded, each term has a phase factor, e _l2T9 ( mi2ai2+m23a23+rra3lQ31 \ and we identify P 12 = § 577712 , P 23 = § 577723 , 
and 531 = §577731 to have 

Ma[ 4>]= yVA[p]e~ i2T i {f,12ai2+p23a23+p3ia3l) . (20) 


Here the weight function Wa [p\ is suppressed exponentially for large p due to the higher power of The fol¬ 
lowing manipulation just goes like the previous case; we can further rewrite as 2 (pi 2 ai 2 + 523023 + p 3 ia 3 i) = 
2 if^f r [ 5 ad-S'adl n bFfi-S'^ d ], where S at \{x) and W^ A are the adjoint counterparts of m and 5 a d is the color charge 

density in the adjoint basis, p a d = p a T a . The symmetric property Ma[( f>aT a } = Ma[U 4> a T a W] for static U requires 
the summation over all {p(.'c)} in a symmetric way here again and we finally arrive at 


M a [</>]= W A [p\ exp 
{p(5D> 



( 21 ) 


III. DENSITY OF STATES 

It is not hard at least to write down a formal expression for the source terms at the present stage. We shall first 
present such an expression and then discuss its physical implication. The source action is deduced from the definition 
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© as 


exp{iSV[A-,p]} = e - i f dix P"' A °- (*) Jvcj)e i S dixpa ' (l, “WM ¥ , A {W cl ’} 

= g -ifd 4 xp a A~(x)jy^ _ 


( 22 ) 


where we defined A f[p\ by 


A f[p] = A4f,a 


Pe 


gfdx^ 


Sp a (x ) 




%] 


(23) 


using that <^„(x) is retrieved by — iS/Sp a (x ) acting on e */rf 4 ® p a <t>a.(x) _ r pj ie phase factor involving A a emerges as a 
result of the variable shift © and, interestingly, this is the only place where A a appears in our expression for the 
source terms. 

This structure of the “weight function” i.e. a function of 8/Sp acting on the delta-functional 8[p\ has been known in 
the color dipole picture i2ql 12711281 ]. If the distribution of p in the present formalism allows only for a quark-antiquark 
pair, the color dipole picture naturally arises. Our expression of N{p\ is not restricted to the color dipole, however, 
and it contains arbitrary distribution of partons also. 

It is interesting that N[p\ has a definite physical meaning as the density of states, namely, the density of allowed 
states under the constraint imposed by a given charge density. That is, roughly speaking, 

-^/■[p] ~ 5 \p( x ) +^( f )] ■ ( 24 ) 

{p(501 


This above expression is, however, not exact because of the nonlinear nature of color charge. Here we shall only make 
sure this interpretation as the density of states in the abelian case in which the above becomes exact. It would be 
helpful to gain some feeling in such a simple case. The counterpart of CD is simply 


q. e 9l dx+ S p‘x) -p g sf dx+ s p S (x) 


6\p] = J2 Wabelian[p] g J d * x ?(*) J dx+ S [ P } , 


(25) 


where Wabeiian [p] allows for p(x) = 0 or p{x) = ±g at each x. In the nonabelian case Jdx + S/5p{x) in the above 
is replace by the nonabelian analogue In W[—i(S/8p a (x))t a ]. Then it is straightforward to reach 1241) by using the 
formula, 


e y d*f(x)=f(x + y). (26) 

Actually the exponential of the p-derivative corresponds to the shift operator caused by the presence of partons. As a 
result of a shift in delta-functional, one parton propagation at x supplies one unit of the contribution to color charge 
density. Hence one can regard as a natural nonabelian extension of the definition of the density of states. 

We shall shortly check the gauge invariance of Sw[A~, p] defined by ( 1221 ) . Under the gauge transformation, 

A~ ^v(^A~ p^VpV\ (27) 

the source term involving the gauge fields is not gauge invariant and acquires an additional inhomogeneous contribu¬ 
tion, 


—iSy Jd A x 2tr[p(x)A (a;)] = —- Jd 4 x 2tr [p(a;){U^(x)9 V(x)}] . 


(28) 


Here we chose to use the color basis in the fundamental representation but the same argument holds in the adjoint 
representation as well. The point is that the density of states is also gauge variant because the color charge itself is 
gauge variant. To investigate the transformation property it is convenient to go back to the first line of and then 
the transformed density of states is read as 


A f[VpV ] ] = 


J V<j>e i ^ x2tT[VpV ^ { - x)] M F ,A\W 4 ‘]. 


(29) 
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Under the variable change 0 by the same V here, tris invariant and thus remains unchanged, while 

the exponential factor in front of it returns to a form close to that before the gauge transformation with discrepancy 
by the inhomogeneous part. We can immediately read the additional term as 


A f[V pV t] = exp 



d A x2ti:[p(x){V\x)d V(a;)}] 


m 


(30) 


which exactly cancels GSI out. Therefore the source terms GU are certainly gauge invariant as a whole, though 
each part is not invariant. The essence is that we do not have to deform A~ to be gauge invariant by means of the 
Wilson line. The inhomogeneous terms resulting from the gauge transformation V involve not A~ but only V and 
p a . Therefore one may as well expect such terms to be eliminated by a factor given only in terms of p“, the role of 
which is actually played by the density of states. In Ref. jTi| the Wess-Zumino term is also characterized as a special 
solution to the transformation property m . It would be an interesting future problem to look for a relation between 
A f\p\ and the Wess-Zumino term, which is not clear at present. 

Although Sw [A ~, p] defined by 1221 ) is certainly gauge invariant and applicable for calculations of the quantum 
fluctuation, if we define the perturbation theory around the classical solution, we do not have to face with nonlocal 
interactions associated with the Wilson line in the source terms as seen in Refs. [DIE!- This would be great advantage 
in a practical sense. We point out that this fact might have been partially observed in the recent calculations of the 
quantum evolution [2£j]. 


IV. MANIFESTLY GAUGE INVARIANT FORM 


It is possible to convert the source terms into a manifestly gauge invariant form written in terms of the Wilson 
line. We started our discussion with the generating functional 0 as a sum of all the contributions with all p a {x). 
We can perform the p-summation not relying on the shift operator formula as argued previously but instead using 
another formula, 





tC = 0 


(31) 


and then we will acquire a different and still equivalent expression for the source terms. After we perform the 
summation over p(x) 7 the source part corresponding to (12211 becomes 


A4f,a 

'p e 9f dx+ jpk^ ta 

J d 4 x p a A a ( x ) 

= A4f,a 

p e ig Jdx +A a ( x)t a 




p=0 



(32) 


Once we recall our results cs and m for the explicit form of the determinant, we can easily deduce from the 
definition 0 the alternative representation of the source terms, 


exp{iSV[A , p]} = Wf[p] exp-| —— Jd 3 x tr[p(ir) In W[A ](x)] 


in the fundamental representation with p = p a t a or 


exp{iS w [A ,p]} = Wa[p] expj--^ yrf 3 a;tr[p ad (f)lnW a d[A ](£)] 


(33) 


(34) 


in the adjoint representation with p ad = p a T a . It should be noted that in this representation the summation over 
p a {x) gives way to the summation over p[x) which derives from the integration over matter B. It has a significant 
meaning that p{x) appearing in (13311 and (1341) should be distinct from p{x). In fact, as we will discuss in Sec.[0 p(x) 
is, so to speak, the classical source if p(x) is referred to as the quantum one. Besides, since p(x) characterizes the 
distribution of fast-moving partons, p(x) is discrete as we have seen in (IT!!!) and 111911 . 

One may have thought that the fundamental and adjoint representations should come about in accord to the quark 
and gluon propagation respectively. This discrepancy is, however, only seeming and not essential at all. We would 
emphasize that the logarithmic form is so peculiar that we should not regard it as merely one variant of the gauge 
invariant generalization nor a simple alternative of ~ tr[pW[A - ]]. As argued also in Ref. [l^, the logarithmic form 
does not depend on the representation, that is, 

-tr[plnIV[A - ]] = Tptr [p ad In Wad [A - ]] • 

9 3 g L 


( 35 ) 











In our calculations we can easily check this from pi2ai2 + P 23 <i 23 + P 31&31 = 3 (p\ai+ P2CI2 + P3CI3) in Sec. Ill Bl From the 
representation independence, it is obvious that the logarithmic form is pure imaginary on its own. This makes a sharp 
contrast to the conventional form ~ tr[p0[.A - ]]. If the source terms are given by ~ tr[pkF[^4 _ ]] in the fundamental 
representation, the complex conjugate ~ tr[pI 0 [A _ ]] would be necessary to make the source terms pure imaginary as 
they should. In the present case of the logarithmic form, on the other hand, the source terms are always imaginary, 
which can be directly seen also from 

jtr[plnIF[A - ]] | = tr[pl lnI0[A - ]] = — tr[plnIF[A - ]] . (36) 

Here we used the hermiticity pi = p and the unitarity 10 [A - ] = IT -1 [A - ]. Consequently the discrepancy between 
EH and 01 resides only in the shape of the weight function Wf,a[p]- 


V. COMMUTATION RELATION 


The color charge operators obey the commutation relation [p a ,p b ] = — igf abc p c , which can be easily confirmed by 
the explicit expressions of p a in terms of quarks and gluons and the canonical quantization conditions. From the 
theoretical point of view of the quantum field theory, [ p a , p b ] = —igf abc p c is necessary for consistency between Gauss’ 
law and gauge invariance, that is, the secondary constraint in the nonabelian gau ge theory (i.e. Gauss’ law) is of the 
first class if we use Dirac’s nomenclature for quantization of singular systems I30l| . 

In the path-integral formalism the operators in the expectation value are time-ordered from the left (large x + ) to 
the right (small a: + ). Therefore the expectation value of the operator commutation relation should be expressed as 
the time-ordered product; 


([P a {x + ), p b (x + )}) = (T(p a (x + )p b (x + - 77 ) - p b {x + )p a {x + - 77 ))) = (p a (x + )p b (x + - 77 ) - p b (x + )p a {x + - 77 )) 


■> 0 + 


7)^0+ 


(37) 

The purpose of this section is to check, in our description given in the path-integral formalism, that the right-hand 
side of the above leads to —igf abc (p c ) as it should. 

In our formulation this noncommutative property follows from the fact that the ^-dependence in the determinant 
appears only through the Wilson line 10 as we have already seen in (ITT1) or llTHl and fHl or (l2ll) . It 
makes no difference whether 10 is given in the fundamental or adjoint representation for the purpose to confirm the 
commutation relation. 

Let us consider the expectation value of color charge density evaluated with the canonical generating functional 
.Zce [p]■ It should be noted that (p a ) can be nonvanishing if it is computed with Zqe[p ], while ( p a ) = 0 after the 
summation over p because of gauge invariance. In this sense, the expectation values (• • •) below are to be considered 
as ones taken at fixed p in the context of the color glass condensate. 

The expectation value of color charge density is easily expressed with the help of the identity derived from the 
integration of a total derivative f <hf){d/d<f> a ){- ■ •) = 0 inside z?ce[p]- It immediately follows 


(P a (*)) 


S\nM[W^} \ 
54>a(x) / 



S In Ad [ 10 ] 




and one more derivative in (f> b from infinitesimally smaller x + yields 


(38) 


(p a (x + )p b (x + - 77 ) - p b {x + )p a {x + - 77 )} ^ o+ = 9 2 (p^^^[W^ x+ (x)[t a ,t b ]W^ + _ T {^].^ 5 i3 \x-y) 

• . £abc / _c / \ r(3) (: 


(39) 


This is exactly what we expect to have. The point is that the subtraction remains nonvanishing only when two 
derivatives act on the same 10 twice which gives rise to the commutation relation between the algebra. 

An independent proof of the commutation relation satisfied with the source terms taking the generic form of (t22|) 
has been provided in Ref. [15]. Let us here discuss another form EH or EH- The charge density p(x) apparently 
would not satisfy the commutation relation. The important point is that, in view of p a {x) can be equivalently 
replaced by iS/SA~(x ) acting on e lSw ^ A 1. In other words, the color charge density operator p a (x) is expressed as 
the derivative in A~(x). (Note that <j> a (x) could be equivalently replaced by — iS/6p a (x ) in Sec. IIVI I In this case it 
is not p{x) but iS/5A~ (x) acting on the source terms which fulfills the commutation property. In this sense p(x) is 
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classical, while iS/6A~{x ) that is p a (x) actually is quantum. BesidesMt should be mentioned that the covariantly 
conserved charge is not p(x) but what is inferred from —SSw/SA~(x) jy|, which also suggests that not p(x) but p{x) 
is the quantum source to be considered. It is almost trivial to check that 


(P a (x)) 


g5W[W^~]] \ 
SAa (x) / 



5S w [W[A-\ 

SW[A-]ij(x) 


\Wt,x+ 


(x)t a W x +_ T 



(40) 


which is the covariantly conserved charge density and one more derivative in A a gives 


{P a {x + )p b {x + -rj) - p b (x + )p a {x + -r])) 


. J 5S w [W[A- 

= l9 \mA- 


r/—>0+ - \dWiA~\ijix) 

=-igf abc ( P c (x )) 6^(x~y). 


[W TtX+ (x) [t a AW x+ My)\u) S {3 \x-y) 


(41) 


From this analysis we can conclude that the noncommutative property has already been incorporated in the framework 
of the JIMWLK equation via the nonlocal source terms written in the Wilson lines even if there is no Wess-Zumino 
term. 


VI. SUMMARY 

We discussed the expression of the source terms which preserves the gauge invariance and the noncommutative 
nature of color charge density operators. It can be written, on one hand, in the form consisting of a naive source term 
with the density of states. Gauge dependence of the naive form is canceled by the transformation property of the 
density of states. Since the interaction vertices associated with the source terms are as simple as in the naive form in 
this representation, the perturbative calculations are expected to go effectively. The work in this direction is now in 
progress. On the other hand, the source terms are to be expressed in a manifestly gauge invariant form in terms of 
the Wilson line. The interesting fact is that the results are not ~ tr[pJV[A - ]] as employed widely but ~ tr [pin W [A - ] 
which had been anticipated from Wong’s equation. We showed how the noncommutative relation of color charge 
arises in our framework, that clearly indicates that the nonlocal interactions in the conventional JIMWLK source 
terms have already handled the noncommutativity. We believe that our new formulation provides a useful hint to 
resolve some subtleties and confusions on the source terms in high energy QCD. 
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